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We consider a phenomenological model of inflation where the inflaton is the phase of a complex 
scalar held $. Planck-suppressed operators of 0{f^/Mpi) modify the geometry of the vev ($) at 
first order in the decay constant /, which adds a first order periodic term to the definition of the 
canonically normalized inflaton (j). This correction to the inflaton induces a fixed number of extra 
oscillatory terms in the potential V ~ 9’’. We derive the same result in a toy scenario where 
the vacuum {3>) is an ellipse with an arbitrarily large eccentricity. These extra oscillations change 
the form of the power spectrum as a function of scale k and provide a possible mechanism for 
differentiating EFT-motivated inflation from models where the angular shift symmetry is a gauge 
symmetry. 


I. INTRODUCTION 

Single-field, slow-roll inflation models that generate de¬ 
tectable levels of primordial gravitational waves require 
a super-Planckian field displacement for the inflaton [1]. 
The constraints on the primordial power spectrum from 
Planck [2-4] remain broadly consistent with the predic¬ 
tions of large classes of these models, particularly those 
with V ^ (jP for 0 < p 2. Knowledge of the de¬ 
tailed predictions of these models in the range p ^ 2 
has become increasingly important as estimates of the 
tensor-to-scalar ratio r have improved [5-7]. 

However, Planck-scale effects should substantially alter 
the simple polynomial relationship V (jP over regions of 
field space that exceed A(() > 1. Monodromy is a mecha¬ 
nism that protects against these dangerous corrections by 
repeatedly wrapping the trajectory of the inflaton around 
a sub-Planckian region of a higher dimensional space, 
with an angular parameter 9. The mechanism relies on a 
corresponding change AV for each winding 9^9 + 27r, 
with the global potential V{(j)) defined by patching to¬ 
gether branch cuts. Large field excursions Acf > 1 are 
then safe because the displacement in the “natural” pa¬ 
rameter space remains small. 

Importantly, monodromy inflation predicts detectable 
oscillatory features in the primordial power spectrum of 
curvature perturbations P(;{k) [8, 9]. If the angular field 
9 is coupled to gauge fields, then non-perturbative in- 
stanton effects modify the potential V ^ V -|-5sin(0i). 
In the absence of monodromy, this oscillatory contribu¬ 
tion can drive natural inflation [10]. Combined with a 
monodromy effect in V, this contribution imprints oscil¬ 
lations in the power spectrum with a frequency that is re¬ 
lated to the periodicity of the underlying theory through 
the axion decay constant /. While the size of these os¬ 
cillations is constrained by the WMAP [11] and Planck 
data [2, 4, 12-14], their detection would provide an un¬ 
precedented glimpse into the monodromy structure of the 
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high energy theory. Therefore, it is important to under¬ 
stand the exact nature of these oscillations and to identify 
any mechanisms that might change this prediction. 

In this paper we study an alternative source of oscilla¬ 
tions in the primordial spectra in monodromy inflation, 
where the angular degree of freedom 9 is the phase of a 
complex scalar $ with a vacuum expectation value (vev) 
of ($) = p. Adding non-renormalizable terms to the ac¬ 
tion for $ alters the shape of the vacuum so that p —>■ p{9) 
and the effective decay constant obtains an angular de¬ 
pendence. These higher order terms only contribute to 
the potential as Unr ~ O (/®), which is much smaller 
than unity for sub-Planckian decay constants / ^ 1. 
However, they change the form of the canonical kinetic 
energy term for <I> in the vacuum at first order in /. Con¬ 
sequently, the canonically normalized inflaton cj), which 
depends on the structure of the kinetic energy term for $, 
obtains a quasi-periodic dependence at where 

m > 5 is the energy dimension of the non-renormalizable 
operator. 

Therefore, the leading order gravitational operators 
with TO = 5 add first order oscillatory corrections to 
any term in the scalar action that depends on the an¬ 
gular coordinate 9. If the shift symmetry in the vacuum 
9 —> 0-|-const, is broken, then the non-renormalizable op¬ 
erators add a set of sinusoidal functions to V that have 
arbitrary relative phases. These new terms are indepen¬ 
dent of the scale of any instanton effects derived from 
gauge couplings and result in a superposition of high 
frequency power spectrum oscillations. These features 
must be taken into account when fitting the predictions 
of monodromy inflation to data. 

Additionally, we consider a case where the structure 
of the vacuum {(j)) is shaped like an ellipse with an arbi¬ 
trarily large eccentricity. The coupling constant for the 
oscillatory terms are suppressed by a factor of at least 1 /3 
compared to the angular symmetry breaking scale, even 
for a highly elliptical vev. A flat prior probability on the 
eccentricity P{^) ~ Z^[0,1], which corresponds to a nat¬ 
urally large correction to the vev, gives sinusoidal terms 
in P^(fc) that are suppressed at the order of 0(10“^) 
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compared to the contribution from the chaotic potential 
V ^ QP. This gives power spectrum features that are still 
safely inside the current bounds from Planck data [12], 

Adding Planck-suppressed terms to the Lagrangian is 
motivated by an effective description of monodromy in¬ 
flation at the energy scale p My>\. However, we do 
require the existence of global symmetries in the low- 
energy theory, which is a non-trivial assumption that 
makes obtaining large-held inflation simple. A more rig¬ 
orous description of the low-energy theory might require 
signihcant hne-tuning in order to obtain a hat potential 
across super-Planckian held ranges. Consequently, we 
assume that an angular symmetry breaking term is tech¬ 
nically natural in a held theory that contains a global 
angular shift symmetry [15] and not include operators 
that would spoil the chaotic inhation behavior of the po¬ 
tential V ^ 9 ^ . We take a phenomenological approach 
and remain agnostic about a specihc high energy scenario 
that realizes this symmetry. 

Monodromy arising from a gauge symmetry is a ubiq¬ 
uitous feature of high energy theories. In particular, 
string theory has many axions, which are angular pa¬ 
rameters that are given a monodromy term by couplings 
to huxes [8, 9, 16-21]. The monodromy potential for 
string axions V ^ (jP can have a variety of exponents, 
with p = {2/3,1,4/3, 2} giving predictions that match 
the Planck posterior probabilities for the ACDM -|- r 
model within approximately l-2cr [21]. While the primor¬ 
dial spectrum resulting from monodromy inflation with a 
broken global shift symmetry has an oscillatory contribu¬ 
tion from each Planck-suppressed operator, the simplest 
models with a gauged axion symmetry give a coherent, 
sinusoidal signal in log k. In principle, these two signals 
can be distinguished via cosmic microwave background 
or large-scale structure data. 

Furthermore, although a large number of angular pa¬ 
rameters 9i could contribute to inflation in principle, 
we only consider the single-field case. While many-field 
models can have simple predictions [22-24], the variance 
of the predicted values of observables for monodromy in¬ 
flation can be undefined for p < 3/4 [25]. Since small 
exponents p are observationally favored, we will leave an 
analysis of many-field monodromy to future work. 

This paper is organized as the following: Sect. II dis¬ 
cusses how non-renormalizable operators alter the shape 
of the vacuum for a complex scalar with a Mexican hat 
type potential. Sect. Ill calculates the canonically nor¬ 
malized inflaton for an arbitrarily elliptical vev. Sect. IV 
estimates the shape of the features in from this 

mechanism both analytically and numerically. Sect. V 
compares the predicted spectra from Planck-corrected, 
EFT-motivated monodromy models to high energy axion 
monodromy scenarios. Sect. VI summarizes the findings. 


II. MODIFYING THE GEOMETRY OF THE 

VACUUM WITH IRRELEVANT OPERATORS 

We start by considering a Lagrangian for a complex 
scalar field, 

$ = (I) 

where the renormalizable terms in the action satisfy the 
canonical Mexican hat shape with 

= ( 2 ) 

We assume that the dimensionless real parameter A is 
0(1) and that the decay constant is / ^ I. The vev 
(‘J’) = / is independent of the phase 9 and is sub- 
Planckian for small decay constants. In the vacuum, the 
radial component of $ has the mass = 4A/^. If 9 is 
coupled to gauge fields, then non-perturbative instanton 
effects induce a sinusoidal correction to Eq. (2), which 
gives rise to natural inflation if / ^ I [10]. 

However, monodromy inflation can be generated in the 
more natural / ^ I limit by explicitly breaking the 9 —>■ 
0-1-const, global symmetry. We do this by adding a term 
to Eq. (2) like 

VAsb=5 0^ (3) 

where the subscript ASH denotes angular symmetry 
breaking. We assume that the form of this term remains 
stable over regions |$| ^ I, which is satisfied in the vac¬ 
uum. The domain of 9 can be analytically extended to 
9 G (— 00 , 00 ) by piecing together each of the branches 
of the Mexican hat, which have a change in potential 
energy of AU = g(27r)^ for every winding. This sym¬ 
metry breaking mechanism has also been invoked as a 
means to dynamically lower the value of the cosmologi¬ 
cal constant [26] and is well motivated from effective field 
theory independently of the monodromy mechanism we 
study here. 

With Eq. (3) large-field inflation can be realized by the 
canonically normalized, pseudo-Nambu-Goldstone boson 
(j) = 9 f when |<I>| = /, if the initial angular displacement 
is sufficiently large 0o ^ lOMpi//. For p = 2 this gives 
an effective inflaton mass of = 2g/f^, which is less 
than the Hubble scale ~ Vasb independently of g 
and is sub-Planckian if g f^. To match the amplitude 
of scalar perturbations as estimated from Planck [7], the 
symmetry breaking parameter should be approximately 
g ^ although the exact value depends on 9o and 

P- 

To ensure that the mass of the radial direction is nrip > 
p[ when the pivot scale leaves the horizon approximately 
A* ~ 55 e-folds before the end of inflation we further 
require 

( 4 ) 
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where we have removed 0(1) terms in the second approx¬ 
imation. A super-Hubble radial mass keeps the evolution 
close to the vacuum ($) = / during inflation, which is 
an important requirement for this paper. ^ 

For a concrete example, Refs [28-30] obtain mon- 
odromy in the phase of a complex scalar with a super¬ 
gravity model where the global U(l) symmetry of the 
Kahler potential is broken by the superpotential in such 
a way that it can give a variety of inflationary potentials 
including Eq. (3). Refs [31-34] have similar helical 
potentials with a range of phenomenological predictions. 
Although realistic scenarios must account for analytical 
control over the symmetry breaking term at high ener¬ 
gies, as well as the stability of a global minimum for 
reheating, we do not worry about these considerations 
here. 

In addition to the Mexican hat shape of Eq (2) we can 
also add any operator that has the form 


— C 



( 5 ) 


0.2 

0.1 

■s 0.0 

Cl. 

- 0.1 

- 0.2 


-0.3 -0.2 -0.1 0.0 0.1 0.2 

p cos 6 



where c is a dimensionless complex coupling, A is the 
cutoff scale, and 

C>^ = ($*’«-h $«). (6) 


FIG. 1. Value of (^) when corrected by Eq. (8) in comparison 
to /. The parameters are A = 1, / = 1/4, and jcj = 1. We 
show the cases (red) g = 1, n = 2; (blue) g = 3, n = 1; and 
(yellow) g = 5, n = 0. 


The operator Om has energy dimension 

m = 2n + q (7) 

and breaks the translational symmetry 9^9 + const, 
in Eq. (2), leaving a residual discrete symmetry 9 —>■ 
9 -I- 27r/g. All operators of this type are generally ex¬ 
pected from effective field theory considerations at energy 
dimension m > 5. With the cutoff scale set to A = Mpi, 
we interpret these operators as gravitational corrections 
to the renormalizable theory at sub-Planckian energies.^ 
The impact of operators like Eq. ( 6 ) has been well stud¬ 
ied in the context of the Peccei-Quinn axion [35], where 
9 is coupled to the QCD gluon field and setting the vac¬ 
uum (9) —>■ 0 dynamically solves the strong-CP prob¬ 
lem. Experimental constraints on CP violation require 
(9) < 10“^° and non-renormalizable terms like Eq. ( 6 ) 
force (9) ^ 0, giving a generally CP-violating vacuum in 
the absence of extreme fine-tuning in the relative phases 
of the operators [36-40]. 


^ Ref. [27] has shown that the dynamics of the heavy radial direc¬ 
tion can reduce the speed of sound of the inflaton fluctuations 
and suppress the amount of gravitational waves. In order to iso¬ 
late the predictions resulting from Planck-suppressed operators 
we do not include these effects in this paper. 

^ While the operators in Eq. (6) are built with Fourier series, we 
could use an alternative basis, e.g., the parametrization <I> = 
X + i'4’ Fas operators Vm = cx”"4’'^ that are also periodic in 9 with 
a similar phenomenology to Eq. (5). However, we ignore those 
bases that have higher order terms that do not leave a residual 
angular symmetry, as well as any possible derivative terms. 


Typically, for inflation driven by the slowly rolling field 
9, an operator Vm >5 is negligible since it is at most an 
0(f^/A) correction to V in the vacuum ($) = = /■ 

However, in the absence of Planck-suppressed operators 
like Eq. (5), the local minimum ($) = pmin is stabilized 
by the |$|^ term in the Mexican hat potential (2). Con¬ 
sequently, the operators O 5 alter the minimum of V in 
the radial direction of 4) at first order in /. Operators 
of higher energy dimension m > 4 change the location of 
the vacuum as 


Pmin 


f 



(y) cos (q9 + 6) 



( 8 ) 


where S is the phase of the coupling c and we have ignored 
terms O {(p — f)^) and terms higher order in /. Fig. 1 
shows pniin for different values of n and q, with the decay 
constant / set to a large value that has been exaggerated 
in order to visualize the change in shape of the vev. 

If we assume that p = Pmin, in addition to the negligi¬ 
ble contributions from Vm> 5 , the deviation of the vacuum 
geometry from / in Eq. ( 8 ) affects the kinetic energy of 
the complex scalar as 





2 


+ P 


2 

min 


d^9d''9. 


(9) 


We define a scalar field ()> with a canonical kinetic energy 
in the usual fashion: 




/ 


d9 



dp 


min 


dO 


2 


■ ■ 
r mm’ 


( 10 ) 
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which is valid up to an arbitrary constant. 

For simplicity we focus on a single leading order non- 
renormalizable operator with energy dimension m = 5. 
The derivative term in Eq. (10) for this operator is 
^ 0{f^) and can be effectively ignored compared 
to the contributions from Up to 0{f^) the inte¬ 

grand reduces to 


de 


1 — /i ^ cos {q9 + 5) 


+ 0(f), (11) 


where we have defined the dimensionless constant 


M = 


5 

8 



( 12 ) 


which we expect to be 0(1) although it could be larger 
if A <C 1. Assuming that the scale of ^ does not disrupt 
the perturbative analysis with respect to /, the inflaton 
can be approximated as 


</>«/ 


9 — sin {q9 + 5) 

9 ^ 



(13) 


which we can invert order-by-order in / to get 


9 


f qA 


sin 





(14) 


In the limit that /x —>■ 0 we recover the standard definition 
0 = </>//. 

Consequently, in the low-energy description we can fix 
P = Pmin E^nd replace all instances of 9 in the scalar sector 
of the action by the relationship given in Eq. (14). Defin¬ 
ing the first and second terms in Eq. (14) as 9 = 9 q + 69i, 
the potential can be expanded like 


Vi9) = Vi9o) + V'{9o)S9i. (15) 


For the combined Mexican hat action (2) and angular 
symmetry breaking potential (3), the effective potential 
at tree level with the Planck-corrected inflaton in Eq. (13) 
becomes 


V = 9 



gppf 

qA 



V 



FIG. 2. The portion corresponding to 0 < 0 < Stt of the 
modified monodromy potential with one additional Planck- 
suppressed operator of energy dimension m = 5, with n = 0 
and g = 5. 


III. A NON-PERTURBATIVE ELLIPTICAL 
VACUUM 

For sub-Planckian field values the corrections to the 
shape of the vacuum from non-renormalizable operators 
like Eq. (6) are naturally small. However, we show in 
this section that the corrections to the potential V can 
be kept perturbatively small even when the vev is altered 
non-perturbatively. 

For a simple example with a non-circular vev we con¬ 
sider the case where (4>) is an ellipse with semi-major axis 
/ and eccentricity which are related by the semi-minor 
axis b as 


e" = i- 



(17) 


plus terms higher order in /. Fig. 2 shows this potential 
for q = 5. Since it is the most conservative choice we use 
a Planckian cutoff scale A = Mpi = I for the rest of the 
paper. 

The second term in Eq. (16) does not come directly 
from the Planck-scale operators (6), but rather from the 
substitution 9{(j)) into Eq. (3), which is due to the modi¬ 
fied shape of the vacuum in Eq. (8). The phase 6 is im¬ 
portant, since any of the possible oscillatory terms from 
Planck-suppressed operators and any oscillations com¬ 
ing from instanton effects generally have different rela¬ 
tive phases, which we discuss in Sect. VB. With p = I 
and /X ~ 0(1), the ratio q/f fixes both the scale and 
frequency of the new oscillatory term compared to the 
angular symmetry breaking scale. 


Instead of the Mexican hat potential (2), we assume that 
U($) is structured so that /3 —>■ /3 -I- const, keeps the po¬ 
tential energy constant and setting ^ t 0 restores the 
9^9 + const, symmetry. We leave the exact form of 
this potential implicit since we will only worry about 
the dynamics when $ = ($). In analogy to the per¬ 
turbative corrections to the vacuum in Sect. II, we keep 
a canonical kinetic energy term for $, which breaks the 
/3-translational symmetry for the Lagrangian when ^ yf 0. 

We parametrize the complex scalar 4) in terms of a 
radial coordinate p and an angular coordinate j3 as 

$ = A cosh (p -|- x/3). (18) 

Surfaces of constant p are ellipses with foci at ±A and 
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FIG. 3. The canonically normalized inflaton as a function of 
the elliptical angular variable /?, where the vev has semi-major 
axis / and eccentricity The dependence is invertible 
but non-elementary. 


we define the vev so that it follows one of these ellipses, 

(^') = Pmin = tanh“^ 0^ . (19) 

We assume that the effective mass in a direction locally 
orthogonal to ($) is sufficiently high so that we can treat 
the system as if it were one dimensional at low energies, 
with p non-dynamical. Near the vacuum we can therefore 
fix the location of the foci at A = f 

The canonical kinetic energy term for <i> in the vacuum 
becomes 

= {f sin^ p + cos^ /3) a^/35'^/3. (20) 

As in Eq. (10) we define the canonical inflaton as 


is composed of a linear and an oscillatory component, 
which Fig. 3 shows are approximately independent even 
as ^ —>■ 1. Motivated by this, we decompose (j){(3) in 
Eq. (22) into a linear and oscillatory portion as 


f 


= 'yi j3 + oscill. 


where Eqs (22) and (23) give the pre-factor 


71 = 

TT 




1 

e-i. ■ 


(24) 


(25) 


This coefficient is bounded in the range 2/7r < 71 < 1 
for 0 < ^ < 1. To determine the form of the oscillatory 
term in Eq. (24) we examine the limiting behavior of the 
elliptic integral. 

For small eccentricities ^ 1 we can expand the ellip¬ 

tical integrals £, which appear in 71 and in Eq. (22), to 
get 


y « (f)sin2/3 + (!l(C"). (26) 

We can compare this directly to the perturbative case 
of Eq. (13) and see that the functional form of these 
expressions matches in this limit, as expected. Similarly, 
for large eccentricities ^ 1, Eq. (22) becomes 

0- |^'0^^sin2/3 + O(C2_l). (27) 

Therefore, Eq. (27) has an identical functional form for 
the oscillatory term as the small-eccentricity expression 
Eq. (26) at first order in Since Fig. 3 shows that 
intermediate eccentricities have similar behavior, a close 
approximation to Eq. (22) for all eccentricities 0 < ^ < 1 
is 


^ — J \lP P + b'^ cos^ (3 dp, 


( 21 ) 


0 

/ 


7i /? + 72 sin 2/?, 


(28) 


which can be expressed in terms of the incomplete elliptic 
integral of the second kind £ as 


P = fVl^£ 


e 1 

^2 - i_ ■ 


( 22 ) 


Fig. 3 shows this relationship for a range of eccentricities 
0 < 5 < 1. Although the function (j){P) is bijective, 
its inverse is not elementary due to the special function 
£. However, we can use the properties of £ to obtain an 
approximation to P{4>) and calculate the potential energy 
of the canonical inflaton cj). 

For /3 > 0, £ is quasi-periodic with period tt, satisfying 


£ y] = 2n£ [tt; y ], (23) 

and £{0;y) = 0 where n is an integer and y is an ar¬ 
bitrary parameter. This relationship ensures that 4>{P) 


where we have defined the coefficient 

72 = Vl-^Af (29) 


in terms of the change in the elliptical integral, which is 


A£ = £ 


'tt ^2 


'tt ^2 

[ 4 ’e-lJ 

2^ 

[ 2 ’e-lJ 


(30) 


The coefficient 72 for the oscillatory term is bounded in 
the range 0 > 72 > —0.21 for eccentricities 0 < ^ < 1. 

Eq. (28) is valid for any ellipse and is not restricted 
to a perturbative expansion in Although this expres¬ 
sion could be inverted numerically to obtain /3(^f)), the 
inverse expression /3(^) is transcendental. An analyti¬ 
cal inversion, as in Eq. (14), can only be made once a 
suitably small parameter is identified. Rather than re¬ 
stricting our analysis to small eccentricities, which would 
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FIG. 4. The inflaton 0 in terms of the angular variable j3 
when the expression (28) is inverted both numerically and 
approximately with Eq. (33). This shows the worst-case sce¬ 
nario with maximum eccentricity ^ > 1, with good agreement 

between the numeric and approximate solutions. 


replicate the results of Sect. II, we instead solve for 
perturbatively in the parameter 

7(0^--. (31) 

7i 

This ratio is naturally small since the oscillatory term in 
Eq. (28) is sub-dominant to the linear even in the large- 
eccentricity limit (27). Importantly, 7 satisfies the exact 
bounds 

0<7<^(v^-l) <^, (32) 

The separation in scales between the linear and oscil¬ 
latory contributions to (j) justifies an order-by-order per¬ 
turbative solution to the inverse function /3((/<). At first 
order this yields 

/3 « j+ 7 sin , (33) 

where the effective decay constant is 

m = ( 34 ) 

which is fixed for a given eccentricity. Eq. (33) is analo¬ 
gous to Eq. (14) in Sect. II. Fig. 4 compares this pertur¬ 
bative inverse function to the exact numerical inverse in 
the maximum eccentricity case ^ > 1. In this limit the 

ratio 7 = — 72/71 maximizes and the approximation is 
close to the numerical values even for this most extreme 
case. Fig. 4 only shows /3 G [ 0 , 7 r], but the periodicity of 
the elliptical integral in Eq. (23) ensures this approxima¬ 
tion is valid for all values of /3. 

Given the relationship (33) we can replicate the anal¬ 
ysis of Sect. II, substituting /3((/) into the low-energy ef¬ 
fective potential. In particular, since we constructed this 


potential so that it has the symmetry /3 —)■ /3-I-const., we 
assume that the angular symmetry breaking mechanism 
acts as 


Fasb =gp^ (35) 


rather than the expression in Eq. (3). Therefore, the 
elliptical vev yields a correction to the potential as 




1 +P7 



(36) 


to first order in 7 . This expression closely matches 
Eq. (16), which arises from modifying the geometry of 
the vev (<I>) = /9„iin(^*) from adding non-renormalizable 
operators to the Lagrangian. However, the frequency of 
the oscillation term in Eq. (36) is not possible with an 
operator of energy dimension m = 5. Therefore the os¬ 
cillatory term above is not Planck-suppressed. Instead, 
large elliptical deviations from a circular vev should come 
from relevant or marginal operator. We leave the exact 
construction of this type of scenario ambiguous. 

This toy model demonstrates that even non- 
perturbatively modified vacuum geometries can still yield 
easily categorizable alterations to the basic shape of the 
potential. As long as we can identify an angular param¬ 
eter { 0 , /?,...} that varies over some finite range in the 
absence of monodromy, then the relationship 4>{0) in the 
modihed vacuum should also be periodic and invertible, 
since Fourier series like Eq. ( 6 ) are the natural basis with 
which to build operators over a compact subset of the 
real numbers. However, one might expect that being able 
to identify perturbative parameters that give a tractable 
perturbative inverse function 9{(f>) is model dependent. 


IV. POWER SPECTRUM FEATURES FROM 
NON-CIRCULAR MONODROMY 


Changing the shape of the vev in Sects H and HI has 
resulted in a modified monodromy potential of the form 


4 "(/>) = 9[j 


p r 


l+p K (-Ism I — 


(37) 


where the constants /, q, and n are related to the parame¬ 
ters in the potentials of Eqs. (16) and (36) in the obvious 
way. The ratio g /sets the amplitude of scalar per¬ 
turbations and K controls the oscillatory contribution to 
the primordial spectra. We estimate the power spectrum 
resulting from Eq. (37) using the 5N formalism. This 
relates the curvature perturbation (/ on equal density hy¬ 
persurfaces to the field perturbations 54> on spatially flat 
hypersurfaces via 


C = N' Scj) (38) 

to first order in Scj), where N is the number of e-folds 
between the flat and equal density hypersurfaces and 
N' = dN/dcj). 
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We approximate the power spectrum whenever a mode 
k crosses the horizon by ignoring any non-standard sub¬ 
horizon evolution for the Fourier transformed field values 
S(j)k- This gives the standard result for the power spec¬ 
trum of scalar perturbations 

'Ps^ik) = ^, (39) 

where Vk = V{(j)k)- From Eq. (38) the primordial curva¬ 
ture power spectrum is 

V^{k) = {N'fVs4,{k). (40) 

Using the slow-roll approximation, the number of e-folds 
of expansion can be calculated as 



k ! kf. 



where 4>c is the value of the field at the end of inflation 
and 4>k is the field value when k = aH. The derivative of 
Eq. (41) is 



P 

+qK 



(42) 


where we have disregarded terms that are 0{k^). The 
third term in the square brackets can be ignored for (j)k 
/, which is the case for those modes that leave the horizon 
within a few e-folds of the pivot scale /c* = 0.002 Mpc“^, 
which has <()* ~ 10 for p ~ 1.^ 

We estimate the horizon crossing value (j>k of the infla- 
ton field in Eq. (42) using only the non-oscillatory portion 
of Eq. (37), i.e., with k —> 0. This is only a rough ap¬ 
proximation to the slow-roll dynamics for (/) and results 
in a value for N' in Eq. (42) that artificially suppresses 
the oscillatory terms, since we are ignoring a component 
to (j)k that is in phase with these terms. A more accu¬ 
rate estimate of (j)k might be obtained by replicating the 
techniques^ of Refs [9, 41] or using generalized slow-roll 
techniques [42-44]. However, this approximation suffices 
to understand the general behavior of the power spec¬ 
trum as a function of scale. 

Setting K = 0 reduces Eq. (41) to 





(43) 


Since the power law term in Eq. (37) gives large-held 
inflation, we ignore any contributions to N that result 


® The derivatives Brjic/dipk vanish for adiabatically evolving scalar 
fields. 

^ This becomes more difficult when the angular symmetry breaking 
term does not have the exponent p = 1. 


FIG. 5. Primordial power spectrum 'Pc(fe) from the poten¬ 
tial (37), normalized by the amplitude of the power spectrum 
at the pivot scale fc,. The parameters are / = 10“^, k = 10“^, 
A* = 55, p — 2/3, and q = 2, which describes an elliptical 
vacuum of eccentricity ^ = 0.28. The solid orange line is 
the 5N approximation and the dashed blue line is the exact 
numerical value. 


from (pc throughout the rest of this section. For the pivot 
scale, Eq. (43) gives 


P, « x/2^. (44) 


We relate pk to the field value when the pivot scale 
leaves the horizon A* « 55 e-folds before the end of in¬ 
flation by 


k _ IkHk _ — Af f 

k:^ y (pifi J 


(45) 


where we have used the slow-roll approximation ~ V/ 
and have related the scale factors at horizon crossing to 
the scale factor at the end of inflation by ak = ace~^. 
Using Eq. (43), this can be expressed as 



P 1 1 

O log — = A* - log 
2 (Pk 


A 

fc* 


(46) 


Modes that leave the horizon within a few e-folds of each 
other have only a small change in cpk- Consequently, we 
solve this equation for (pk while ignoring log(pif/(pk in or¬ 
der to understand the behavior of P(^{k) only near fc*. 
The solution to Eq. (46) in this limit is 



(47) 


This expression is valid only for k/k^ ^ e^* and when 
(pk depends only weakly on the oscillatory terms in V. 

Substituting Eqs (42) and (47) into Eq. (40) gives an 
estimate of the primordial power spectrum. Fig. 5 com¬ 
pares this approximation to an exact numerical calcu¬ 
lation for the potential (37) using the publicly available 









inflation solver MultiModeCode [45-48]. MultiMod- 
eCode numerically solves the Klein-Gordon and mode 
equations for (/) on sub-horizon scales and is exact to first 
order in 5(j)k- The chosen parameters for Fig. 5 corre¬ 
spond to an elliptical vev with eccentricity ^ = 0.28, 
which results in sinusoidal oscillations in Vc^ as a function 
of logk/k^. The analytic solution is a close approxima¬ 
tion, although the amplitude of the oscillations is sup¬ 
pressed relative to the exact calculation, as expected. 
For smaller amplitude oscillations with |k| < 10“^ the 
approximation becomes a better fit. 

Constraints on primordial tensor modes from Planck 
and BICEP2 [2-4, 6 ] restrict p ^ 2, although this is 
model dependent. To match recent studies that constrain 
oscillations in the primordial spectrum [11, 13, 14, 49], we 
generally expect |k| < 10 “^ and / < 10 “^ for the fiducial 
case of p = 1. When Eq. (37) is generated by Planck- 
suppressed operators, the upper limit on |k| corresponds 
to p < 0.5 with p defined as in Eq. (12).^ 


V. DIFFERENTIATING PLANCK-CORRECTED 
AND GAUGED OSCILLATIONS 

A. The toy case 

In Sects II and III we considered an effective descrip¬ 
tion of monodromy inflation where a global angular sym¬ 
metry 9^9 + const, is explicitly broken and irrelevant, 
Planck-suppressed operators. This differs from many ax- 
ion monodromy scenarios that come from explicit high 
energy considerations, such as string axions [51] or ex¬ 
tranatural inflation [52], where the translational symme¬ 
try in 6 * is a gauge symmetry. 

Typically, the field 9 in high energy theories is identi¬ 
fied with an axion, which is obtained by integrating gauge 
fields over cycles in compact extra dimensions. The resid¬ 
ual gauge symmetry for this scalar severely restricts the 
allowed structure of the Lagrangian and precludes opera¬ 
tors like those in Eq. ( 6 ) from perturbatively altering the 
shape of the vacuum ($) once the gauge is fixed. How¬ 
ever, the gauge symmetry in these cases requires that 4) 
is coupled to gauge fields, which induces an oscillatory 
term to the potential via non-perturbative instanton ef¬ 
fects, which take the form 


hyp = 5 np cos 



(48) 


^ The data can allow much wider ranges for the parameters, since 
the posterior probabilities for the model parameters can be highly 
non-trivial. Furthermore, large wavelength oscillations in V(^(k) 
can be confused for models that have large (positive or nega¬ 
tive) running [50], due to the finite range of multipole moments 
that can be measured at high accuracy in the cosmic microwave 
background. 


We treat this as the prototypical potential for gauged 
monodromy inflation, although the model might be dif¬ 
ferent in more specific cases, particularly those that 
incorporate a natural variability in the oscillation fre¬ 
quency [31, 32, 34, 41]. 

In principle, signals in the primordial spectrum due 
to sinusoidal terms like Eq. (37) should be distinguish¬ 
able from the features derived exclusively from non- 
perturbative instanton effects like in Eq. (48). The dif¬ 
ference is most obvious for the case where the vacuum 
geometry is modified by a single Planck-suppressed term 
in the Lagrangian and the potential is closely approxi¬ 
mated by Eq. (16). Most obviously, the coupling con¬ 
stant for the sinusoidal term varies as a function of (j) for 
p 7 ^ 1. For p < 1 modes k that leave the horizon closer 
to the end of inflation feel a relatively larger effect from 
the oscillatory portion of Eq. (16) than those that leave 
the horizon earlier.® The resulting signal in the power 
spectrum for these smaller modes may be detectable in 
the correlation functions of large scale structure in large- 
volume surveys [53-56], which are able to constrain the 
amplitude of small modes /c ^ fc*. 

The non-renormalizable oscillations also have a larger 
frequency than the instanton terms by an integer multi¬ 
ple. 


UJq = 


q 

f 


vs 


wnp 


1 

7 ’ 


(49) 


where q G {1,3,5} for the leading Planck-scale correc¬ 
tion. This gives features in P(;{k) that are of higher fre¬ 
quency than features from Eq. (48) for hxed values of 

/■ 

Furthermore, if we assume that all dimensionless cou¬ 
plings in Eq. (16) are exactly unity, then the coupling 
for the new sinusoidal term from V 5 as compared to the 
leading order polynomial term is 


P k- f ^ Z. 
q [sj ujg' 


(50) 


With the angular symmetry breaking exponent p fixed, 
this only has the frequency uiq as a free parameter. The 
relative amplitude of oscillatory features in the primor¬ 
dial spectrum P({k) that result from the instanton and 
non-circular vev effects is then directly tied to their rela¬ 
tive frequencies. However, the amplitude of power spec¬ 
trum oscillations that result from Eq. (37) is independent 
of those coming from Eq. (48). 


B. Power spectrum with multiple oscillatory 
components 

In general, the analysis in Sect. V A is overly simplistic, 
since we should allow for all possible non-renormalizable 


For 0 <IC / and p < 1, the sinusoidal term dominates over the 
diverging pre-factor (0//)^~^ so that V —>■ 0 here. 





9 



We again assume that the total contribution of the oscil¬ 
latory terms above only negligibly affects the relationship 
4>k{(t>*) in Eff- (47). Therefore, the power spectrum P({k) 
can be estimated by simply substituting N' from Eq. (52) 
into Eq. (42). The power spectrum then roughly corre¬ 
sponds to a power-law with a superposition of sinusoidal 
terms, each of which has a form as in Fig. 5. Fig. 6 
displays the power spectrum for the superimposed signal 
with two different choices of couplings /ii and phases 6i. 
This relatively complex primordial signal differs substan¬ 
tially from the easily distinguishable oscillations expected 
from the potential (48) alone. 


k / 


VI. SUMMARY 


FIG. 6. Power spectrum when all lowest-order Planck- 
suppressed operators of 0{f^) are included. The angular sym¬ 
metry breaking exponent is p = 1/3 and the decay constant is 
/ = 0.1. The relative phases between the different oscillatory 
contributions and their overall couplings are chosen indepen¬ 
dently from the distributions Si ~ W[0,27r] and pi ~ W[0,1]. 
Two different selections of fii and Si are plotted. 


operators to contribute to the modification of the vac¬ 
uum, including derivative terms that we have ignored 
throughout this paper. However, adding only those 
Planck-suppressed operators of energy dimension m = 5 
that are described in Eq. ( 6 ) changes the relationship 
9{(j)) in Eq. (14) to 

« = 7 + + ( 51 ) 

where the sum is over the three values of qi S {1,3,5} 
that are allowed from Eq. (7). The coefficients are de¬ 
fined as in Eq. (12) for each of the qi, where the phases 
6i are assumed to be uniformly distributed. Eq. (51) is a 
perturbative inverse function that is valid to 0{p) and 
requires that the sum of the individual couplings is sub¬ 
dominant to the linear term 9 = (j^/f- This marginally 
lowers the range of decay constants for which this expres¬ 
sion is valid in comparison to Eq. (14). Potentials with 
multiple oscillatory terms have been studied in the ab¬ 
sence of monodromy in Refs [57-60] and with monodromy 
in Refs [61, 62]. However, the tight relationships between 
the coupling constants and frequencies in Eq. (51), com¬ 
bined with the limited number of Planck-suppressed op¬ 
erators of energy dimension m = 5, severely reduce the 
effective dimensionality of parameter space for this sce¬ 
nario in comparison to previous work. 

The general expression (51) modifies the potential (37) 
so that it contains a summation of out-of-phase oscilla¬ 
tion terms. This similarly affects the derivative of the 
number of e-folds as 



In Sect. H we studied the low-energy effective descrip¬ 
tion of monodromy inflation resulting from the evolution 
of the angular component 0 of a complex scalar field $. 
We started with the only renormalizable potential that 
gives a non-zero vev ($) = / and assume / <C Mp\. 
Motivated by effective field theory we included higher 
order terms with energy dimension m > 5. Although 
these terms are suppressed at order we showed 

that they modify the geometry of the vev so that it has 
an angular dependence that is affected more severely: 
(4>)^/[l + O(/( 0 ))]. 

When we restrict $ = ($), the new geometry for the 
vacuum adds an angular dependence to the canonical ki¬ 
netic energy term for $. In the absence of irrelevant 
operators, the canonically normalized scalar field that 
drives slow-roll inflation is </) = f 9, but acquires a sinu¬ 
soidal dependence when these terms are included. This 
correction is first order in / and the inverse function is 
given schematically by 0 ^ (j)/ f + 0{f) sin(u;(/)). 

In Sect. HI we looked at more radical modifications to 
the shape of the vacuum. We parametrized the vev as an 
ellipse with arbitrary eccentricity 0 < ^ < 1 and found 
that the elliptical angular coordinate variable /3 is related 
perturbatively to the inflaton (jj in terms of the small pa¬ 
rameter 0 < 7 (^) < 1/3, which is true for all possible 
eccentricities. This parameter describes the ratio of lin¬ 
ear and oscillatory contributions to the elliptic integral 
£{n-,y), which is guaranteed to be small since y{^) < 0 
for any ellipse. This gives a description of the effective 
potential in the vacuum that is analogous to Sect. H, 
without requiring a perturbatively small modification to 
the circular Mexican hat vev. 

Since the relationship between the angular parameter 0 
and the inflaton cj) has a Hrst order correction, a term V ^ 
9P that breaks the symmetry 9 —>■ 9 +2tt is also corrected 
at first order in /. The modified potential obtains one or 
more sinusoidal terms, which each contribute oscillations 
to the primordial power spectrum P(;{k). In Sect. IV 
we calculated the power spectrum resulting from these 
types of corrections and found that they are oscillations 
in log k. 

Without a reason to exclude Planck-suppressed opera- 
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tors from the theory, the existence of multiple oscillatory 
terms in the potential is completely general. As discussed 
in Sect. V, the P({k) from many such terms should be 
closely approximated by the linear superposition of sig¬ 
nals with varying phases. However, the individual signals 
have a close relationship between their frequencies, which 
is constrained to a few possible integer values due to the 
energy dimension of the non-renormalizable operators. 
Furthermore, the overall coupling constant for the terms 
that determine the size of oscillations in P(;{k) can be 
related to the oscillation frequency, if the dimensionless 
parameters in the theory are 0{1). 

The EFT-motivated description of the low-energy the¬ 
ory of a complex scalar contrasts with those cases where 
monodromy inflation results from the simplest non- 
perturbative breaking of a gauge symmetry. In principle, 
the effects of the modified vacuum mechanism presented 
here would distinguish a high energy axion monodromy 
theory from monodromy arising in an effective low-energy 
description. Consequently, searches for non-sinusoidal 
features in the primordial power spectrum might even¬ 
tually allow us to extract a Planckian signal from the 
primordial spectrum. Conversely, the absence of compli¬ 


cated oscillations in F’<;(fc) would constrain the allowed 
form of Planck-suppressed operators in monodromy in¬ 
flation, adding to existing constraints on high dimension 
operators from local non-Gaussianity [63]. However, ob¬ 
taining significant evidence in favor of a non-sinusoidal 
feature in the primordial spectrum remains a challenging 
problem. 
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